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Abstract 

We examine the effect of a more realistic Earth matter density model which takes into account 
of the local density variations along the baseline of a possible 2100 km very long baseline neutrino 
oscillation experiment. Its influence to the measurement of CP violation is investigated and a com- 
parison with the commonly used global density models made. Significant differences are found in the 
comparison of the results of the different density models. 
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1 Introduction 



The first results from KamLAND collaboration have recently been announced. The data demonstrate 
the disappearance of the reactor v e at a high level of confidence [Q and hence corroborate the oscillation 
solution to the solar neutrino problem. Furthermore, the measurement excludes all but the Large Mixing 
Angle oscillation solutions (LMA) [21 EI - I n the fit of all solar neutrino data, prior the KamLAND result, 
including the neutral and charge currents, elastic scattering of both day and night data, the sign of 
the Am^ has been determined to be positive better than 3cr for the solar mixing angle, # S oiar < tt/4 
0J|S]. Among all the theoretical and phenomenological implications 0, the LMA solution establishes a 
favorable condition for a determination of the leptonic CP violation (CPV) phase as a further probe 
of new physics in long baseline neutrino oscillation experiments [S]. 

Because of the smallness of the mixing angle #13 HJ, the signal of CPV effect will not be large in 
general. This makes the measurement of the CP phase a challenging task. Hence a detailed estimate of 
the various possible theoretical uncertainties and experimental errors will be crucial for the extraction of 
the CP phase. In particular, the matter effect has to be properly delineated in very long baseline (VLBL) 
oscillation experiments. The presently available Earth matter densities are modeled globally averaged 
values. They are given as functions that depend on the Earth radius only. However, Earth density is not 
a spherically symmetric function independent of the longitudinal and latitudinal coordinates. There are 
local density variations and can have abrupt density changes from place to place, radially or at different 
longitude and latitude. Although the average density model may be suitable for a variety of purposes, one 
wonders if it is adequate for an accurate extraction of the lepton CP phase by VLBL experiments. We 
can identify two specific questions which may affect the outcome of VLBL experiments and to which we 
have to look for answers: One question is how do we assign a realistic error to the model matter density? 
The other question is how do we estimate the effect of local matter density deviations from the available 
average? Unless we find satisfactory answers to these questions, we can not be sure that the errors in the 
extraction of the CP phase from VLBL is under control so that we can assign a good confidence level to 
the value of the CP phase obtained. 

We have addressed the first question in a recent publication JD] where we proposed a set of density 
profiles which are randomly distributed around the average density to simulate the way that Earth matter 
density is determined. This approach provides a way to estimate the error, induced by the uncertainty 
of Earth matter density, on the CP phase determination. Other approaches have subsequently been 
proposed and a summary of several approaches can be found in [TT] . 

In this paper we address the second question. We consider a more realistic matter density function 
along a specific baseline so that we can use a concrete example to examine the question. We will again 
focus on the recently approved high intensity proton synchrotron facility near Tokyo, Japan, i.e., the J- 
PARC (Japan Proton Accelerator Research Complex) ^21- We consider a VLBL, with and beams 
from J-PARC to a detector located near Beijing, China, as depicted in Fig. A preliminary study of 
the possibility of such a VLBL experiment, which we called H2B, can be found in Refs. J3| and |14|. 
In ^2 and subsequent studies a number of physics issues have been investigated: physics potentials 
of H2B, relevant backgrounds and errors JSJj the effects of Earth matter density uncertainties '1,0: , 
and the feasibility of measuring CP violation and atmospheric neutrino mass ordering in two joint LBL 
experiments JB]- Other studies on the matter effect on CP can be found in |17| . 

In Section 2 we discuss briefly the Earth density models and propose an alternative, more realistic 
density model for H2B. Section III presents a quick review of the approach of Ref. JU| f° r dealing with 
uncertainties of Earth matter density. Section IV discusses briefly the general statistical and other errors 
in CP measurements in LBL. A brief discussion is presented in Sec. V. 

2 Earth density models and matter density distributions 

In looking for the Earth matter effect, we are usually provided with some global model of Earth matter 
density, e.g., the PREM ^H] or AK135 [T§|. All presently available Earth density functions are not 
directly measured but obtained using a limited set of geophysics data which are analyzed by means of 
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an inversion procedure. The density so obtained is a function of the depth from the Earth surface and 
any longitudinal and latitudinal variations are ignored. Consequently, in a given density model, the same 
density profile will be given for all baselines of the same length irrespective of their locality. All oscillation 
experiments of the same baselines length will have the same matter effect. Furthermore, since the matter 
density is a symmetric function along the baseline the matter effect will be the same when the neutrino 
beam source and detector sites are interchanged. Because of these simplified features, the existing density 
models are inherently of limited level of precision for VLBL. To improve the precision we have to know 
the specific density profile for a given VLBL. Or we have to establish the fact that the effect of the density 
variation is within the tolerance of the uncertainties that exist for the experiment. 

To demonstrate the effect of local density variations we make a detailed examination of the mass 
density profile along the baseline of H2B. Figure [21 shows the Earthquake P-wave velocity perturbation 
around the AKf 35. The color codes the size of the deviations from AKf 35. By mapping out the deviations 
along the baseline, we can obtain a more realistic density profile for H2B. We shall call this density profile 
for H2B the H2B-FTL density function. 2 We will describe below how H2B-FTL is obtained. Three density 
profiles along the H2B baseline are shown in Fig. El PREM by the dotted curve, AKT35 by the dashed 
curve, and H2B-FTL together with its typical error bars by the solid curve. The three density profiles 
provide a significant range of density variations which allow us to investigate the effect of Earth density 
on the determination of neutrino oscillation parameters in VLBL, in particular the CP phase. 

Similar to PREM and AKf 35, H2B-FTL is based on the traveling time of earthquake waves. However, 
it is concentrated on the local region of the H2B baseline and therefore involves a much larger set of 
available data along the baseline than either PREM or ALf35. Moreover, since it is a 3-dimensional 
density model, it contains significant information in the longitudinal and latitudinal directions, while 
PREM and AKf35 are f-dimcnsion density models which provide information only along the radial 
direction. Hence H2B-FTL can better represent the actual density along the H2B path than either 
PREM or AKf35. 

The density function of H2B-FTL is related to that of AK135 via the following relation |2()l 121 j . 

1 Sv 

PH2B-FTL = (1 + ~, JPAK135) (1J 

K V 

where the AK135 density function, PAK135, can be found in |19| . v is the P-wave velocity and Sv is the 
P-wave velocity correction to AK135 |22|- The geophysics consideration of the H2B path gives k = 0.2. 
The ratio of the P-wave velocity correction to the P-wave velocity is the sole geophysics input in the 
corrections to AKf 35 and is given in terms of a large set of discrete data on the various positions along 
the baseline [201 - As shown in the scale at the bottom of Fig. [21 Sv/v varies from +2.3% to -8.3%. 

As shown in Fig. [31 in several regions, H2B-FTL deviates from AKf 35 (PREM) beyond the usually 
cited allowed variation of ±2% (±5%). The strongest deviations can be recognized as the red and blue 
regions of Fig. [21 The red and blue regions represent respectively negative and positive corrections to 
AKf 35. As Fig. ^ displays, in its path from J-PARC to Beijing, a neutrino will go through the upper 
Earth mantle which exhibits plastic properties. It will first experience the Japan island crust (blue), 
the asthenosphere raised by Pacific slab in collision with the European- Asian slab (red), the normal 
asthenosphere under the Sea of Japan (red), the bottom of the West-CoSon-Man of west Korea peninsula 
(red), and the Bo Hai Sea of China (red). As shown in Fig. [21 global density models, such as PREM and 
AKf 35, can have significant deviations from the actual mass density profile of the H2B baseline. 

The actual variance of the H2B-FTL is a complicate function along the baseline. For simplicity and 
to be conservative, we will ignore its position dependence and just take the maximal square root variance 
|20j to represent the density variance, 

Sv 

a( — ) = 0.003 (2) 

V 

Then we have 

O"(,0H2B-FTl) PAK135 1 , Sv 

= T a \ J 

PH2B-F TL PH2B-FTL k V 

2 The H2B-FTL density file is based on the work of one of the authors and his geophysics research group 1201 . 
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= Q.Q15 x PAK135 

PH2B-FTL 



(3) 



which lies between 0.75% and 3%. It should be mentioned that the H2B-FTL density profile consists 
of a huge data sample, hence the discretization size of the baseline in the path integral can be as small 
as 40 km, while the discretization size of PREM and AK135 along this baseline is 200 km. A smaller 
discretization size can reduce the error which is another factor that contributes to the higher level of 
precision of H2B-FTL. 



3 Errors in CP violation measurements 

We review briefly the method of ^U] an( A define our notation. We are interested in quantifying the 
possible error caused by Earth density uncertainties. As usual, we define a CP-odd difference of neutrino 
oscillation probability functions, 

D(S C v,N e (x)) = P a p(5 CP ,N e (x))-P & p(6 CP ,N c (x)), (4) 

where Sq p is the CP phase, N c (x) the electron density distribution function along the baseline x, and P a p 
{Paf}) the oscillation probability of v a — » Vp (y a — > Vp). The electron density function N e {x) is related 
to the Earth matter density p(x) by the Avogadro number, N\, and the electron fraction, Y c , through 
the usual relationship, N c (x) = NaY c p(x). The dependence on the neutrino energy, mixture angles, and 
neutrino masses is suppressed. 

We estimate the possible error from density uncertainty by following the formalism of |10| . 



SD(S GP ,N e (x)) ee yj < {{P af} - P & p)- < (P a0 - P~) >) 2 > 

= \f(8(P a p)f + (5(P a p))\ (5) 

where < • • • > denotes a weighted average of a matter density dependent quantity. It is defined as a 
path integral of the quantity, along a given baseline, over an ensemble of possible variations of Earth 
matter density profiles, N c (x), weighted by a logarithmic normal distribution functional F[N e (x)] for 
non-negative quantities |1U| . e.g., 

<P a p> = J[VN c (x)]F[N c (x)]P a0 (S CP ,N c (x)). (6) 

Therefore we interpret the average matter density N c (x), such as PREM and AK135, as a weighted 
average over all samples of possible density profiles, 



N c (x) ee < N c (x) >= j [VN c (x)]F[N c (x)}N c (x). (7) 
The matter density uncertainty is given as usual by 

(8) 



a(x) ee ^<N?(x) > -(N (xj)' 



The level of precision of a density model with a given N(x) and cr(x) can be measured by 

r(x) ee 9 
N e (x) 

To do the functional integral, the integration path along the neutrino baseline is discretized according 
to the available geophysical information. For the H2B baseline, the discretization size for H2B-FTL is 
40 km, which is much smaller than the 200 km discretization size suitable for PREM or AK135. The 
smaller discretization size is helpful in reducing the errors in the functional integration. 
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We remark that an alternative parameterization of the Earth matter density uncertainty is to give the 
average matter density a fixed deviation, i.e., taking the density function to be N'(x) — (1 ± r'(x))N(x), 
where r'(x) is the density uncertainty. Conventionally N(x) is given by PREM or AK135 and r'(x) is 
respectively 0.05 or 0.02. As discussed in [10) . this parameterization will lead to a larger uncertainty in 
the extraction of the CP phase than the present approach. 

The above formulation can be readily adopted to the measurable quantity of the event number. Unless 
noted otherwise in the calculation of event number, we assume 6500 interacting muon neutrinos for H2B. 

For the numerical calculation, we adopted the following values for the solar and atmospherical neutrino 
mass square differences and the corresponding mixing angles: 3 

Am soiar = 5.0xlO" 5 eV 2 , tan 2 6>i 2 = 0.42, 

Am 2 tm = 3.0xlO" 3 eV 2 , sin 2 2(9 23 = 0.99. (10) 

Since the CPV effect is proportional to 6*13 , a larger value of sin 2 26*i3 will be more favorable for its 
measurement. The CHOOZ bound is sin 2 2013 < 0.15 [5]. In Fig. 0] we plot the CPV event number vs 
the CP phase. We show several different values of sin 2 26>i3 and the beam energies. One can see from the 
upper panel Fig. 01 that, for this number of testing muons there will be little sensitivity in the CP phase 
measurement in H2B if sin 2 26*i3 < 0.01. However with a higher number of testing muons the sensitivity 
can be increased. In the following calculation we will use sin 2 26*13 = 0.08 for illustration. However, for 
comparison, we will also show some results for sin 2 26*i3 = 0.01. 

To select the appropriate neutrino energy, we have examined the oscillation probabilities of 1.5, 4, 4.5, 
5, and 8 GeV. As shown in the lower panel of Fig. 4.5 GeV is the optimal energy which will be used in 
all subsequent calculations. It should be noted that Fig. 0] employs the density distribution H2B-FTL, 
but different density model does not change the optimal energy and the sensitivity in 6*13. 

To examine the difference among the different density models we define the following quantity which 
provides a gross measurement of the "pure" CP effect, 

AD(5 CP ) = D(5 CP ,N c {x)) - D(6 C p = 0, N c (x)). (11) 

Then we can define a quantitative measure of the significance of the CPV signal by the ratio of the "pure" 
CPV effect and the corresponding quantity which also contains the matter effect 4 , 

_ AD(S C p) f . 

Odcn(OGP) = JJ^TT „ -, y. (12) 

dD{d C p, N e (x)) 

A larger Sden(<5cp) gives a stronger CP signal relative to the matter effect, we require that Sden(^cp) be 
significantly greater than 1. In Fig. [SJ we plot <Sd e n(<$cp) against the CP phase Sep- Indeed H2B-FTL 
gives a signal better than AK135 and PREM for Sep away from and tt where the CPV effect vanishes. 
To compare the H2B-FTL with other density model, e.g., the AK135, we define 

5D(5 C p) = D(5 CP ,N e (x)) 

(H2B-FTL) (AK135)- 

(13) 

We can define a significance measure of the CPV effect with the difference of density models, 

- AD(Scp) , . 

Sden(dCP) = -~— -• (14) 

dD(dcp) 

We plot <Sdon(<5c) in Fig. We see that the difference in density models is always larger than the CP 
difference. This again underlines the fact that a realistic density is necessary. 

In Fig. 0we plot the "pure" CP effect, AD{§cp) given in Eq. fTTJ) . as a function of the CP phase 
for the three density models together with the error bars for H2B-FTL. In general, AK135 is 3<r away 
and PREM is 6<r away from H2B-FTL for a given CP phase, indicating significant differences between 
the different density models. 

3 The recent KamLAND best fit Q] gives Am^ tm = 6.9 X 10~ 5 and 0.86 < sin 2 26>i 2 < 1.0. The value of sin 2 20i 2 used 
in the present work corresponds to the lower limit of the KamLAND value. 

4 In Ref. 1101 (see Eq. (17) there) we used the inverse of S ( jen('5cp) as defined in Eq. 1121 . Since AD(Scp) vanishes in 
the absence of CPV, we find the present definition is more convenient to graph. 
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4 Statistical and other Errors in CPV measurements 



In this section we present the result of our study of the statistical errors, the background effects, and 
the systematic uncertainties, and contrast them with the effect of density uncertainties. Our treatment 
of statistical errors, background effects, and systematic uncertainties follows that of Ref. ^5]. They are 
represented by their respective error factors denoted as / (statistical), r (background), and g (systematic). 
Denoting the variance of their combined effect by ctsbs j we can define the significance of measure of the 
CPV effect with respect to this combined variance, 

S'sbs = _^2e2 (15) 

0"SBS 

Obviously Ssbs is only meaningful for the case of a sufficiently large number of interacting v^s. For the 
numerical calculation we take the statistical factor / = 0.02, the background factor r — 0.01, and the 
systematic factor g — as used in ^j- We show in the upper left corner of Fig. a representative error 
bar for this combined error. 

In Fig. Elwe plot both S S bs, Eq. HSJl, and S dcn , Eq. for sin 2 20 13 = 0.08 and 0.01 and for two 

total number of interacting v^'s, 6500 and 650. We see that for several hundred interacting i/^'s, no CPV 
effect is expected to be measurable. For several thousands of muons, we have a chance to observe the 
CPV effect. With the CPV variable we consider here Earth matter uncertainty can be ignored safely if 
H2B-FTL model is employed, but the significance will be worse if other density model (Sden) is adopted. 

5 Discussions 

In this paper we investigated the potential error caused by the uncertainties of earth matter density. We 
also consider the other errors, such as the statistical error, in the context of intended H2B physics. We 
found that the more realistic density model, H2B-FTL, the matter density variation induces a rather small 
error which allows a meaningful separation of H2B-FTL from AKI35 and PREM. With CPV variable we 
used, i.e., the CP difference, the statistical and other errors generally dominate over the difference of the 
density models. Therefore, a different approach of the CP measurement, such as that used in jTBj has to 
be investigated. 
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Figure 1: A schematic diagram of the H2B baseline which consists of a neutrino superbeam from J-PARC 
which is located about 60 km northeast of Tokyo, Japan to a detector near Beijing, China. The longitude 
and latitude of the neutrino source and target are indicated. The baseline is 2100 km. 
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Figure 2: A 2-D plot of the earthquake P-wave data based on which the H2B-FTL density model is 
constructed. The scale from Tokyo to Beijing is the longitude in degrees and the scale perpendicular to it 
is the depth measured from Earth surface in units of km. The neutrino trajectory of H2B is a symmetric 
arc from Tokyo to Beijing, with its deepest penetration of Earth about 90 Km near 128° longitude. 
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1000 
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Figure 3: The vertical axis is the Earth effective matter density profile, where Y e p(x), Y c — 0.494 is the 
electron fraction and p{x) is the matter density in gram/cm 3 . The horizontal axis is the baseline from 
J-PARC (0 km) to Beijing (2100 km) in km. The PREM and AK135 are world average densities and 
H2B-FTL is constructed specifically for H2B. The error bars on the solid curve are the uncertainties 
of H2B-FTL defined in Eq. Note that the densities have been scaled by a factor of the nuclear 

composition, Y c = Z/A ~ 0.494, which is suitable in the region of H2B. 
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Figure 4: The v and D c number difference, corresponding to Eq. |@J with 6500 interacting z/ M , vs the 
leptonic CP phase. The curves in the upper panel are for different values of sin 2 20\^ and those in the 
lower panel are for different energies. The horizontal axis which is the CP phase Sqp is in units of 2n 
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Figure 5: The CPV signal significance S dcn , Eq. O, vs the CP phase. The solid line is for H2B-FTL 
with a 3% uncertainty, the short-dashed line is for AK135, and the dashed line for PREM. The upper 
panel has sin 2 26*i3 == 0.01 and the lower panel sin 2 2#i3 = 0.08. The scale of the horizontal axis is in 
units of 2-7T. 
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Figure 6: The CPV signal significance 5dcn, Eq. 1)14(1. vs the CP phase. The solid curves is for sin 2 20i3 = 
0.01 and dashed curve for sin 2 2#i3 = 0.08. The horizontal axis is in units of 2ir. 
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Figure 7: The CPV difference SD(5 C p), Eq. ©, vs the CP phase for sin 2 20 13 = 0.08. The squares 
are for AK135, the circles for PREM, and the stars for H2B-FTL. The error bars associated with the 
stars are for H1B-FTL. A typical size of the statistical plus background errors is shown in the upper left 
corner. See the relevant discussion in Sec. 4. 



15 




Figure 8: The CPV signal significance vs the CP phase. The solid line is for Sden, Eq. 114(1. The dashed 
lines are for the signal significance Sqbs, Eq. I|15fl . The thin dashed line has 650 interacting and the 
heavy dashed line has 6500 interacting v^. The upper panel is for sin 2 2#i3 = 0.01 while the lower panel 
for sin 2 26»i3 = 0.08 
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